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Abstract. We studied in detail the binary cluster can- 
didates SL538 and NGC 2006 in the Large Magellanic 
Cloud (LMC). This cluster pair is located in the north- 
western part of the large OB association LH 77 in super- 
giant shell LMC 4. A third star cluster, KMHK 1019, is lo- 
cated within 5' from the cluster pair. Our study is the first 
age determination of SL 538 and NGC 2006 that is based 
on CMDs. We derived an age of 18 ± 2 Myr for SL 538, 
22.5±2.5 Myr for NGC 2006, and 16 Myr for KMHK 1019. 
Thus the three clusters are (nearly) coeval. We identified 
Be star candidates and find the same ratio N(Be)/N(B) 
for the components of the binary cluster (12%) while the 
amount of Be stars detected in KMHK 1019 (5%) and in 
the surrounding field (2%) is considerably lower. Since Be 
stars are usually rapid rotators this may indicate intrin- 
sically higher rotational velocities in the components of 
the cluster pair. Also the IMF derived from the CMDs 
shows the same slope for both SL 538 and NGC 2006 and 
is consistent with a Salpeter alMF. An estimation of the 
cluster masses based on the IMF slopes showed that both 
clusters have similar total masses. These findings support 
joint, near-simultaneous formation of the cluster pair in 
the same giant molecular cloud. 

Key words: Hertzsprung-Russel (HR) diagram - stars: 
emission line, Be - stars: luminosity function, mass func- 
tion - star clusters: individual: NGC 2006 and SL 538 - 
Magellanic Clouds 
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1. Introduction 

The existence of gravitationally bound pairs of star clus- 
ters is important for the understanding of formation and 
evolution of star clusters. Since the probability of tidal 
capture of one cluster by another one is very small (Bhatia 
et al. 1991), we can assume that the components of a true 
binary star cluster have a common origin. Star clusters 
form in giant molecular clouds (GMCs) (e.g. Elmegreen 
& Elmegreen 1983), but the details of cluster formation 
are not yet well understood. If the components of a clus- 
ter pair or multiple cluster formed simultaneously or se- 
quentially in the same fragmenting GMC they should have 
similar properties like age, metallicity and stellar content. 

In the Milky Way only a few binary clusters are known, 
though Lynga & Wramdemark (1984) suggest the com- 
mon origin of a group of six Galactic open clusters. Later 
Pavloskaya & Filippova (1989), and more recently Subra- 
maniam et al. (1995), propose the existence of more pos- 
sible Galactic binary clusters and cluster complexes. 

The apparent lack of binary clusters in our own Galaxy 
may be explained in different ways. Subramaniam et al. 
(1995) argue that since we are looking at the Galaxy from 
inside, double clusters may be harder to detect than in 
the distant Magellanic Clouds, where binary clusters can 
easily be detected due to the closeness of their projected 
positions on the sky. The distance to the Galactic clus- 
ters must also be taken into account, but only approxi- 
mately 400 of 1400 open clusters have known distances 
(Lynga 1987). Subramaniam et al. (1995) found 16 Galac- 
tic binary cluster candidates on the base of the Lynga 
catalogue, which corresponds to 8 % of the investigated 
number of clusters. From this they conclude that binary 
clusters in the Milky Way may not be uncommon. 

The evolution of a gravitationally bound pair of star 
clusters depends on the interaction between the compo- 
nents as well as on the tidal forces of the parent galaxy. If 
the tidal field is strong, the binary system will not survive 
for long but soon will get disrupted. From some prelim- 
inary considerations Innanen et al. (1972) conclude that 
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due to stronger tidal forces in the Milky Way a binary clus- 
ter will execute only a fraction of a single orbit around the 
barycentre before its components are detached, but it will 
survive for several orbits in the less dense, less massive 
Magellanic Clouds. Surdin (1991) came to the same con- 
clusion, especially for massive clusters. The investigation 
of binary clusters may help to evaluate the tidal field of 
the parent galaxy. 

Fujimoto & Kumai (1997) suggest that globular and 
populous star clusters form through strong collisions be- 
tween massive gas clouds in high-velocity random mo- 
tion. Shear and momentum of oblique cloud-cloud colli- 
sions lead to break-up into compressed sub-clouds revolv- 
ing around each other, which may form binary or multiple 
clusters. Binary star clusters are expected to form more 
easily in galaxies like the Magellanic Clouds with high- 
velocity random gas motions, whereas in the Milky Way 
such large-scale high- velocity random motions are lacking. 

Bhatia & Hatzidimitriou (1988), Hatzidimitriou & 
Bhatia (1990), and Bhatia et al. (1991), have surveyed the 
Magellanic Clouds in order to catalogue the binary clus- 
ter candidates. The maximum projected centre-to-centre 
separation of the components of a pair was chosen to be 
18 pc, which corresponds to ~ 1.'3 in the LMC. A bi- 
nary cluster with larger separation may become detached 
by the external tidal forces while shorter separations may 
lead to mergers (Sugimoto & Makino 1989, and Bhatia 
1990). In these studies 69 pairs in the LMC and 9 pairs 
in the SMC were identified. Two clusters may appear to 
be a binary cluster due to chance line-up while in fact 
being at different distances within the Magellanic Clouds 
and not gravitationally bound to each other. The number 
of chance-pairs of objects uniformly distributed in space 
can be estimated with a formula presented by Page (1975). 
Taking into account also a non-uniform distribution of star 
clusters (at least for the LMC), Bhatia & Hatzidimitriou 
(1988) and Hatzidimitriou & Bhatia (1990) found that 
statistically 31 pairs in the LMC and 3 pairs in the SMC 
could be explained due to mere chance line-up. As consid- 
erably more pairs have been found, this strongly suggests 
that at least a certain amount of them must be true binary 
clusters. 

While it is difficult to measure true distances between 
apparent binary clusters an analysis of their age and stel- 
lar content can give clues to a possible common origin. 

The star cluster pair NGC 2006 (also known as SL 537) 
and SL 538 is located in the northwestern part of the OB 
association LH 77 in supergiant shell LMC 4, and has a 
projected centre-to-centre separation of 54'.'67 correspond- 
ing to 13.3 pc. This double cluster has already been the 
subject of investigations concerning its binarity. Bhatia 
(1992) and Bica et al. (1996) found from integrated pho- 
tometry that both clusters have the same age. Kontizas et 
al. (1993) analyzed the stellar content and the cores of the 
components using low resolution objective prism spectra 
and integrated IUE spectra. They suggested this cluster 



pair constitutes a true binary cluster, which moreover may 
merge in some 10 7 years. This raises another question, 
which has also been discussed in Bhatia & MacGillivray 
(1988): could mergers of former binary star clusters be re- 
sponsible for at least some of the blue populous clusters 
in the LMC? 

We investigated the double cluster NGC 2006 and 
SL 538 in an attempt to find further affirmation - or disaf- 
firmation - of the binarity of the two clusters. We analyze 
the star density in the clusters and the surrounding field 
(Section 3). For the first time we derive ages for these 
clusters from isochrone fits to colour-magnitude diagrams 
(CMDs) (Section 4), which is a much more reliable age 
determination than using integrated photometry. In Sec- 
tion 5 we investigate the content of Be stars in the clusters 
as well as in the surrounding field. In Section 6 we give a 
summary and conclusions. 

2. The data and the data reduction 
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Fig. 1. F-image of SL 538 and NGC 2006, north is up and east 
to the left. The field of view of our image is 5'8 x 5^8. To the 
south-east of the cluster pair a third star cluster is located: 
KMHK 1019. The location of the Be star candidates (see Sect. 
5) is marked with circles 



An image of SL 538 and NGC 2006 is shown in Fig. |. To 
the south-east of the cluster pair a third star cluster is 
located: KMHK 1019. 

The data have been obtained on 15 December 1992 
with EFOSC 2 at the ESO/MPI 2.2 m telescope at La 
Silla. A 1024 x 1024 coated Thomson THX31156 chip 
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(ESO #19) was used with a pixel scale of (X'34 result- 
ing in a field of view of 5. '8 x 5. '8. All data were obtained 
with the standard Bessell B, V, R, Gunn i and Ha filters 
used at the 2.2 m telescope (see Table [l] for an observing 
log). 

After standard image reduction with MIDAS, profile 
fitting photometry was carried out with DAOPHOT II 
(Stetson 1991) running under MIDAS. 

The photometry was transformed using the Landolt 
standard fields around PG 0218+029, Rubin 149, and 
SA98 (Landolt 1992) observed in the same night. 

We applied the following transformation relations: 

b - B = z B + a B ■ X + c B ■ (B - V) 
v - V = z v + a v - X + c v ■ (B - V) 
r - R = z R + a R ■ X + c R ■ (V - R) 
i - I = zi + aj ■ X + cj • (V - I), 

where X is the mean airmass during observation, capi- 
tal letters represent standard magnitudes and colours, and 
lower-case letters denote instrumental magnitudes after 
normalizing to an exposure time of 1 sec. The resulting 
colour terms Cj , zero points Zj and atmospheric extinction 
coefficients are: 

z B = 2.593 ± 0.049 mag 
a B = 0.140 ± 0.039 mag 
c B = -0.257 ±0.007 
z v = 1.271 ± 0.038 mag 
a v = 0.135 ±0.030 mag 
c y = -0.063 ± 0.006 
z R = 1.170 ±0.050 mag 
a R = 0.081 ± 0.039 mag 
c R = 0.007 ±0.015 
zi = 2.504 ± 0.040 mag 
a/ = 0.073 ±0.032 mag 
a = 0.047 ±0.005 

3. Stellar density in and around the clusters 

In order to investigate the stellar surface density distribu- 
tion in and around the clusters we subdivided the entire 
area into 52 square cells with 20 pixels (corresponding to 
6'.'8) length each. We counted the number of stars found 
within each cell. In Fig. ^ (left) one cell corresponds to one 
pixel of the star density plot. The counting was carried 
out on the ALLSTAR-output tables, which contain the 
coordinates of each measured star. We see an enhanced 
star density between the cluster pair. The apparent stel- 
lar bridge between the components contains as many as 
25% of the stars counted in each of the clusters. To make 



Table 1. Observing log 



Object 


Filter 


Exp. time 
[sec] 


Seeing 
["] 


SL538/NGC2006 


B 


360 


1.5 


SL 538/NGC 2006 


B 


90 


1.5 


SL 538/NGC 2006 


V 


150 


1.3 


SL 538/NGC 2006 


V 


20 


1.1 


SL 538/NGC 2006 


R 


150 


1.3 


SL 538/NGC 2006 


R 


30 


1.4 


SL 538/NGC 2006 


Gunn i 


180 


1.3 


SL 538/NGC 2006 


Gunn i 


20 


1.3 


SL 538/NGC 2006 


Ha 


900 


1.2 



density structures and thus the stellar bridge between the 
two clusters better visible we applied a 3 x 3 average filter 
for image smoothing (see Fig.^, right). Since we expect 
especially high values in the stellar bridge compared to 
the surrounding field, a median filter - which does not 
consider very high or very low pixel values - seems not 
to be appropriate for our purpose. However, we also tried 
a median filter and found no differences concerning the 
features of the resulting images. 

No signs of an increased star density reaching from 
KMHK 1019 towards the binary cluster candidate can be 
seen. 

We investigated the statistical probability of seeing 
such a stellar bridge due to chance density fluctuations. 
For this purpose we removed all stars outside a 60 pix- 
els radius around the components of the cluster pair and 
carried out 50 artificial star experiments with ADDSTAR 
running under DAOPHOT II. Each time we added the 
same number of stars at random coordinates as had been 
removed initially. Star density plots were created in the 
same way as described above. We found 12 plots out of 50 
showing larger star density somewhere between - and con- 
necting - the star clusters. In most cases the stellar bridge 
connecting the components is very thin. Only three plots 
showed stellar bridges wider than one cell. In Fig. || three 
examples of our artificial star density plots are shown: no 
bridge, small bridge, and wide bridge. Our artificial star 
experiments show that in one out of four cases a stellar 
bridge may occur due to statistical density fluctuations in 
the field. Thus, the probability of a stellar bridge which 
is not due to statistical star density fluctuations does not 
reach a confidence level of 95 %. 

4. Deriving ages for the star clusters 

We derived ages of the clusters SL538, NGC 2006, 
KMHK 1019, and the surrounding field star populations 
by comparing our CMDs with isochrones. The isochrones 
we used are based on the stellar models of the Geneva 
group (Schaerer et al. 1993). 
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Fig. 2. Left: Investigating the stellar surface density distribution we see enhanced stellar density between the two components 
while no enhancement between the binary cluster candidate and KMHK 1019 can be seen. We plotted 8 grey values, reaching 
from black (no star per cell) to white (7 stars per cell). Right: To make density structures and thus the stellar bridge between 
the two clusters better visible we applied a 3 x 3 average filter for image smoothing. Note that one cannot use this image for 
star counting as the pixel values are averaged to repress noise 
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Fig. 3. Star density plots based on artificial star experiments. Left: no bridge occurred between the star clusters, middle: a 
small bridge connects the components of the cluster pair, right: the stellar bridge is more pronounced. This last case occurred 
only three times out of 50 star experiments and shows that there is a low likelihood of a pronounced stellar bridge produced by 
random fluctuations (6%) 



4-1- Age determination 

We derived colour-magnitude diagrams (CMDs) of each 
cluster by cutting out a circular area with a radius of 
100 pixels, corresponding to 34" or 8.3 pc, centered on 
the optical centre of each cluster. To derive the CMDs of 
KMHK 1019, which is the smallest of the three clusters, 
we adopted a smaller radius of 80 pixels, corresponding 



to 2772 or 6.6 pc. An estimation by eye based on the star 
density plot (Fig. |J) suggests that no or almost no cluster 
stars are outside this area. 

All CMDs are plotted in Figs. | to @. Each cluster 
CMD has a wide blue main sequence and contains very few 
supergiants. The width of the main sequence is caused in 
part by photometric errors, crowding (seeing 1'.'3) and the 
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Fig. 4. Colour-magnitude diagram of the star cluster SL 538. We adopt an age of 18 ± 2 Myr. The derived age is supported by 
the best-fitting isochrones in all diagrams. Be star candidates (Sect. 5) are marked with crosses 
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Fig. 5. Same as above, but for the cluster NGC 2006. The adopted age is 22.5 ± 2.5 Myr, which is confirmed by the isochrone 
fits in all colours 



presence of Be stars (see Sect. 5). The scarcity of super- 
giants is well within the expected fluctuations for compact 
clusters with few stars. 

Overplotted on the CMDs are the best fitting 
isochrones. We fitted the isochrones such that the su- 
pergiants rather match the blue loops than the quickly 
traversed subgiant branch. A distance modulus of 18.5 
mag (Westerlund 1997) was adopted. The metallicity of 
the young field population of the LMC was found to be 
[Fe/H] ~ —0.3 ± 0.2 dex by various authors (Russell & 
Bessell 1989, Luck & Lambert 1992, Russell & Dopita 
1992, Thevenin & Jasniewicz 1992). We therefore adopted 



Geneva isochrones with Z=0.008 which corresponds to 
[Fe/H] ~ -0.3 dex. 

Galactic field stars contaminate our observed area. 
Ratnatunga & Bahcall (1985) estimate the number of fore- 
ground stars towards the LMC, and in Table || we present 
their counts scaled to our total field of view (5'8 x 5'8) . 

SL 538: We see two red and one blue supergiant in the 
V, B-V CMD of SL 538. The 16 Myr (solid line) isochrone 
fits all three supergiants well but also the 20 Myr (dotted 
line) isochrone fits the red supergiants very well. We adopt 
an age of 18±2 Myr. The same age is found from isochrone 
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Fig. 6. CMD of KMHK 1019, the smallest of the three star clusters and with the lowest number of stars. The resulting age for 
the cluster from the best fitting isochrone is 16 Myr 




B-V [mag] V-I [mag] 

Fig. 7. CMD of the surrounding field. This diagram comprises the mixture of ages of the various field populations: the blue 
main sequence and the supergiants represent the younger populations, the intermediate age populations show up through the 
RGB and the pronounced clump. Note the widening of the main sequence due to Be stars (see Sect. 5). Stars with V > 14.5 
are plotted as smaller dots to keep the isochrones recognizable 



fits to the V,V-I CMD (Fig. |). All isochrones are based 
on a reddening of Eb-v = 0-05 mag. 

NGC 2006: Four supergiants are located in this cluster, 
covering a colour range from B — V = to 1.5 mag. Both 
CMDs (Fig. H) are fit quite well by isochrones with ages of 
20 Myr (solid line) and 25 Myr (dotted line). We adopt an 
age of 22.5 Myr (±2.5 Myr) and a reddening of Eb-v = 
0.05 mag. 

KMHK 1019: This cluster by far is the smallest one 
with the lowest number of stars. Few data points are lo- 



cated in the red clump and the red giant branch (RGB) 
region, and it is very likely that these stars belong to an 
intermediate-age field star population, while the super- 
giants, on which our age determination mainly relies, are 
located in the cluster centre and thus we assume that they 
belong to the star cluster. The main sequence is sparse, 
especially in the upper part brighter than V ~ 16 mag. 
The best fitting isochrones in both CMDs result in an age 
of 16 Myr (Fig. |). The reddening of E B -v = 0.05 mag is 
the same for all fits. 
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Table 2. Number of foreground stars towards the LMC calcu- 
lated from the data of Ratnatunga & Bahcall (1985), scaled to 
our field of view of 5. '8 x 5. '8 



apparent visual magnitude range 



colour range 


13-15 


15-17 


17-19 


19-21 


21-23 


(B — V) < 0.8 


1.3 


2.9 


3.1 


6.4 


6.0 


0.8 <(B-V)< 1.3 


0.5 


2.6 


5.4 


4.3 


7.4 


1.3 < (B — V) 


0.1 


0.6 


3.7 


13.1 


29.3 



The surrounding field: The field population comprises 
a mixture of ages. Apart from a blue main sequence and 
blue and red supergiants, which represent the young field 
populations, the intermediate- age field population of the 
LMC shows up through red giants and the pronounced 
red horizontal branch clump. We are not able to distin- 
guish between distinct young populations, but the plotted 
isochrones represent ages which are supported by corre- 
sponding supergiants. 

The brightest blue supergiants and some of the bright- 
est red supergiants are represented by the 16 Myr 
isochrone (solid line). Also the 25 Myr (dotted) and 32 
Myr (dashed) isochrones are supported by bright blue, 
yellow and red supergiants. 

Several supergiants are traced by the 80 Myr isochrone. 
Note that the redder main sequence stars at V ~ 16 mag 
are candidate Be stars (see Sect. 5). One could easily mis- 
take them for more evolved stars marking an additional 
field population with an age of approximately 100 Myr. 
The stellar density seems to be lower between 80 Myr and 
200 Myr which indicates a possible decrease in the field 
star formation rate. 

Along the 200 Myr isochrone and below, the star den- 
sity is increased along the subgiant branch, again indicat- 
ing enhanced star formation. 

Eb-v = 0.02 mag is a lower limit to the reddening 
of the field star populations and corresponds to the blue 
envelope of the main sequence(s). 

The youngest field population is part of LH77. Our 
derived age of approximately 16 Myr is in good agreement 
with the findings of Braun et al. (1997). 

4-2. Comparison to earlier photometry 

Our study is the first age determination of SL 538 and 
NGC 2006 that is based on CMDs. Previous studies de- 
rived ages based on surface photometry, using different 
aperture sizes. Age determinations based on integrated 
colours are less precise than age determinations based on 
CMDs: Geisler et al. (1997) investigated the influence of 
a few bright stars on the integrated light of intermediate 
age star clusters. They conclude that fluctuations in the 
number of bright main sequence stars and red giants lead 
to shifts in the integrated colours, which affect the age 
determination. In Table we present a comparison to the 



integrated photometry of Bica et al. (1996), Bhatia (1992) 
and our results. 

Though several authors using surface photometry state 
that integrated colours and thus the derived ages are 
largely independent from the aperture radius, Bica et al. 
(1996) and Bhatia (1992) derived quite different ages for 
SL 538 and NGC 2006: Bica et al. (1996) used an aperture 
size of 50" and found NGC 2006 to be the older compo- 
nent of the cluster pair. This is in agreement with our re- 
sults. In contrast, Bhatia (1992) used an aperture radius of 
33" and found SL 538 to be slightly older than NGC 2006. 
Our CMDs exclude ages as young as suggested by Bhatia's 
(1992) aperture photometry. Bica et al.'s (1996) young age 
of only 0-10 Myr for SL 538 again is clearly excluded by 
our data, while the wide age range of 10-30 Myr for NGC 
2006 includes our result of 22.5 ± 2.5 Myr. Note that Bica 
et al.'s (1996) apertures are so large that stars belong- 
ing to the neighbouring cluster are also included in their 
measurement. 



Table 3. Comparison of earlier age determination based on 
integrated colours and ours based on CMDs 



reference 


aperture size 


cluster 


age 


Bica et al. (1996) 


50" 


SL 538 
NGC 2006 


0-10 Myr 
10-30 Myr 


Bhatia (1992) 


33" 


SL 538 
NGC 2006 


12.6 Myr 
7.9 Myr 


this work 


34", CMD 


SL 538 
NGC 2006 


18 Myr 
22.5 Myr 


5. Be stars in 


the clusters 


and the 


surrounding 



field 



Be stars are non-supergiant B stars with variable Balmer 
emission and infrared excess originating in circumstcllar 
disks. Rapid rotation and red/infrared excess of Be stars 
lead to redder colours and widened main sequences (see 
Grebel et al. 1996). The observed brightness and colour of 
a rotating (Be) star depend on its rotational velocity and 
inclination to the line of sight (rotational displacement 
fan, Collins & Smith 1985). 

Using the R—Ha index to detect stars bright in Ha (R 
serves as continuum filter) and B — V as a temperature 
index we can identify Be star candidates (Fig. |J). This 
method was first described in Grebel et al. (1992, 1993). 
In Fig. U a pronounced clump of data points at R—Ha ~ 
mag and B — V ~ mag can be seen. These data points 
represent blue main sequence stars and blue supergiants 
without Ha emission. To the red side of this clump, scat- 
tered around R — Ha ~ 0.1 mag and extended over the 
whole B — V colour range, stars belonging to the RGB and 
red supergiants are visible in the two-colour diagram. Red 
giants and red supergiants can have Ha emission, and in- 
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Fig. 8. Two-colour diagram to detect Be star candidates. 
R serves as a continuum filter to detect stars bright in Ha 
(R — Ha > 0.2), while B — V serves as a temperature index 
to separate the blue stars from red giants and supergiants. For 
more information see text 



deed some data points are located at higher R — Ha values 
(up to ~ 0.35 mag). Our selection criteria for the candi- 
date Be stars are B — V < 0.2 mag and R — Ha > 0.2 
mag. 

The widening of the main sequence of the CMDs due 
to Be stars can be clearly seen in Figs. || to [?], where we 
marked the Be star candidates with crosses. The effect is 
most pronounced in V — I since the I filter has the highest 
sensitivity to the infrared excess of Be stars. 

In Fig. [I] we marked these stars with circles. Our Be 
star candidates are concentrated at the location of the 
components of the cluster pair whereas at the location of 
KMHK 1019 only two such stars are present. 

The impression that the Be star candidates are domi- 
nantly present in SL 538 and NGC 2006 is confirmed when 
considering the ratio of Be stars to B stars. Since we do 
not have spectral classifications we simply considered the 
ratio of all B to Be stars within a magnitude interval of 
V = 14.2 to 19.1 mag. These magnitudes correspond to 
the mean visual magnitudes, at LMC distance, of B III 
to B9V main sequence stars (Zorec & Briot 1991, Table 
3). We find the following values: 



SL538: 
NGC 2006: 
KMHK 1019: 
field: 



N(Be)/N(B) 
N(Be)/N(B) 
N{Be)/N{B) 
N(Be)/N(B) 



0.123; 

0.120 
0.053 
0.019 



+0.145 
-0.078 
+0.233 



+0.022 



The errors are corresponding to 3<t Gaussian errors and 
are calculated using the confidence limits for small num- 
ber statistics from Gehrels (1986). The components of the 
double cluster show the same fraction of Be stars, and the 



sixfold amount found in the surrounding field. The differ- 
ence between the ratios N(Be)/N(B) of SL538 and the 
field is 0.104 which is more than a 3cr-cffcct. The differ- 
ence between SL 538 and KMHK 1019 is 0.070 which is less 
than a 2<7-effect according to the upper confidence limit 
for KMHK 1019. Thus, the Be star content of KMHK 1019 
may be comparable to the cluster pair. 

6. Summary and Conclusion 

Investigating the stellar density around the clusters (see 
Fig. ||, right) we see no signs of increased star density 
reaching from KMHK 1019 towards the cluster pair, but 
we see an enhanced star density between SL 538 and 
NGC 2006. This may indicate a stellar bridge and thus 
gravitational interaction between SL538 and NGC 2006. 
The stars in the bridge are main sequence stars and thus 
they may either belong to the cluster "system" or to the 
young field star populations. However, artificial star exper- 
iments showed that a stellar bridge may also be explained 
(24% probability for a "small" but only 6% probability 
for a more pronounced bridge) by statistical fluctuations 
in the field star density. 

Fitting Geneva isochrones (Schaerer et al. 1993) to the 
CMDs we find the following ages: SL538: 18 ± 2 Myr, 
NGC 2006: 22.5 ±2.5 Myr, KMHK 1019: 16 Myr, youngest 
field: 16 Myr. The three clusters might have formed se- 
quentially as part of the same GMC that formed LH 77. 

Be stars are concentrated in SL538 and NGC 2006, 
and both clusters show the same ratio of N(Be)/N(B). 
This is in agreement with Kontizas et al. (1993) who in- 
vestigated integrated IUE spectra and the distribution of 
spectral types of stars and found that both clusters have 
similar stellar content. Since Be stars are usually rapid ro- 
tators this may indicate intrinsically higher rotational ve- 
locities in the components of the cluster pair. The amount 
of Be stars detected in the surrounding field is consider- 
ably lower. 

An investigation of the IMF of the binary cluster can- 
didate showed that the IMF slopes agree with each other 
within the errors and are compatible with the Salpeter 
value (r = —1.35, Salpeter 1955). Our results for the 
IMF slopes are F = -1.22 ± 0.31 for SL538 and T = 
-1.27 ± 0.32 for NGC 2006. With these values we esti- 
mate the upper limits for the total cluster masses to be 
2300±1100M© (SL 538) and 2300±1200M Q (NGC 2006). 
The similarity of the cluster masses is in agreement with 
the findings from Kontizas et al. (1993). Let us assume 
that the cluster pair indeed is a binary system: Follow- 
ing Kepler's third law and assuming 4600M© for the total 
mass of the binary system and 13.3 pc for the distance 
between the clusters we get an orbit period of « 47 Myr. 
This would mean that the clusters so far have moved less 
than half of an orbit since their formation. 

The similarity of properties (ages, Be star content, 
slope of the IMF and masses) indicates possible joint for- 
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mation and suggests small spatial separation. Our results 
suggest that SL 538 and NGC 2006 are a true binary clus- 
ter, but a final proof is still missing. Radial velocity mea- 
surements would help to understand the kinematics of the 
clusters and may prove or disprove its possible binary na- 
ture, but such data are not yet available. If the stellar 
bridge is real this may imply that the merger process has 
begun. However, our artificial star experiments showed 
that we cannot state whether this bridge is real or not. 
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